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When using Einstein’s equations, there exist a number of techniques for embedding non¬ 
linear structures in cosmological backgrounds. These include Swiss cheese models, in 
which spherically symmetric vacua are patched onto Friedmann solutions, and lattice 
models, in which weak-field regions are joined together directly. In this talk we will 
consider how these methods work in f{R) theories of gravity. We will show that their 
existence places constraints on the large-scale expansion of the universe, and that it may 
not always be possible to consider the Friedmann solutions and weak-field solutions of a 
theory independently from each other. 


1. Introduction 

There are a number of methods that can be used to consistently include non-linear 
structures in cosmological models. The most well-known of these is probably the 
“Swiss cheese” approach, in which spherically symmetric solutions are joined to 
homogeneous and isotropic Robertson-Walker geometrie^ (see Fig. [^a)). An 
alternative approach is to take regions of perturbed Minkowski space, and join 
them together at reflection symmetric boundaries!^ (see Fig. I^b)). The former of 
these approaches has the benefit of admitting exact solutions. The latter benefits 
from admitting space-times with reduced symmetry, and from removing the regions 
of locally homogeneous and isotropic space-time entirely. 



(a) (b) 


Fig. 1. Two methods of embedding non-linear structures into a cosmological model, (a) A Swiss 
cheese model, (b) A lattice model. The solid filled region indicates a homogeneous fluid. 


In the context of Einstein’s theory, the Swiss cheese and lattice models can be 
used to elucidate the links between the gravitational fields of non-linear structures, 
and those of the cosmological models within which they reside. They show that 
the Friedmann solutions of Einstein’s equations are entirely compatible with the 
Coulomb-like gravitational fields that are expected on small scales. One does not 
restrict the other. In this paper we will consider the corresponding situation in the 
/(i?) theories of gravity. We will reach quite different conclusions in this case. 
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2. The f{R) Theories of Gravity 


The f{R) theories of gravity are a generalisation of Einstein’s theory of gravity. 
They are defined by replacing the Ricci scalar, R, in the Einstein-Hilbert action, 
with some function, f{R). Variation with respect to the metric then gives the 
following set of field equations: 


T ^ 

r — _ 

— p ' £ 
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2 (/ - RfR)g^.. + 


( 1 ) 


where subscript R denotes a derivative with respect to that quantity, and where 
is the energy-momentum tensor of the matter fields (defined in the usual way). 
These equations are fourth order in derivatives of the metric, and if one want to 
join together two space-times then the following conditions® are required for their 
union to be a solution of Eq. Q: 

fR[dyK:,t=0, [dyKt=Q. [ 9 . 7 : 6 ]+= 0 , ( 2 ) 

Irr [dyR]t = 0, = 0, (3) 


where dy is the normal derivative to the boundary, "fab is its induced metric, and K 
and are its trace and trace-free parts of the extrinsic curvature, respectively. 
The brackets [...]^ denote the difference in a quantity when evaluated on either side 
of the boundary, and we have assumed here that there are no surface layers {i.e. 
there are no J-function matter fields on the boundary). The first three of these 
equations are the same as required by Einstein’s theory (up to a factor of Jr), while 
the last two are new. These latter equations constitute additional requirements that 
a geometry must obey, in order to be a solution of the field equations. 

In much of the analysis reported here, we will us a post-Newtonian appoximation 
to describe weak gravitational fields in the presence on non-linear structures. The 
leading-order part of the perturbed metric is then given by 

ds^ ~ —(1 — 2(j))dt^ + (1 + 2ij})5ijdx^dx^ , (4) 


where (j) and <C 1 can be interpreted as gravitational potentials. The Taylor series 
approximation to theories that can be expanded around i? = 0 is then given by 

/(i?) = /(O) + fR{0)R + \fRR{m^ + 0{R^). (5) 

If we take /(O) ~ /^(O) ~ IrrIO) 1, then the gravitational potentials become® 

^ ^ + ^JrrR^ and tlj = ^ ^JrrR^ , (6) 

where U, U and R are functions that satisfy = —47rp — /(0)/2, V^C/ = 

— 47rp -I- /(0)/4, and SJrrS/^R — JrR = —Sirp + 2/(0). The reader will note that 
the first two of these potentials are exactly the same as the potentials that result 
from solving Einstein’s equations with a cosmological constant (where /(O) takes 
the place of the usual symbol A). The Ricci scalar, R, plays the role of an additional 
Yukawa potential. 
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3. Swiss Cheese in f{R) Gravity 

Let us first consider Swiss cheese constructions, in theories that can be expanded 
as in Eq. ([^. If we model the region around a point-like mass using the post- 
Newtonian approximation, discussed above, then an application of the junction 
conditions in Eq. (|^ tells us that the radius of the vacuole that surrounds it must 
obeyl^ 

^__ k 

E2 “ /flS S2 + E ^ ’ 

S U,r\i: fRRR,rh /o', 

S - /flE ^ 2/j^ E ’ ^ ^ 

where the boundary is at coordinate distance r = T, from the centre of the vacuole, 
and where fc is a constant. The first two terms in each of these equations behave 
in the same way as the contributions from dust and A in the Friedmann solutions 
of Einstein’s equations. The second term in Eq. Q behaves in the same way as a 
spatial curvature term, and the last terms in Eqs. and ([^ have no counterparts 
in the solutions of Einstein’s equations. 

In we now impose the remaining junction conditions, from Eq. ([^, then we see 
that R r = 0 on the boundary, such that the last terms in each of Eqs. Q and 
([^ must vanish. Recall that these are the only two terms that show any difference 
from the solutions of Einstein’s equations, and one is led to the conclusion that 
the only Swiss cheese models that can be constructed within this class of theories 
are ones that expand at leading order in a way that is identical to the Friedmann 
solutions of Einstein’s equations with A. This is a strong restriction, imposed on 
the large-scale expansion simply by demanding that the theory have a well-defined 
post-Newtonian limit in the vicinity of point-like masses. 

4. Lattice Models in f{R) Gravity 


The lattice models, while mathematical! and conceptually quite different from 
the Swiss chees models, display a remarkable similarity in the final result of their 
construction. In this case it can be shown, as a result of both the field equations 
and the junction conditions, that the position of the boundary must satisfyl^ 






JrZ 


(9) 


where here the boundary is taken to be at coordinate distance z = Z from the 
centre of the lattice cell. This is again for theories that can be expanded as in Eq. 
([^. The right-hand side of the equation above takes the same form as the source 
terms in the Friedmann solutions of Einstein’s equation in the presence of dust and 
A. We therefore have the same result as before: The only large-scale expansion 
that is possible, in the leading-order behaviour of these constructions, is identical 
to that of the Friedmann solutions of Einstein’s equations with A. 
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5. Non-Analytic f{R) Models 

The results derived so far seem quite limiting, but one might question whether 
they are a result of assuming the theory can be expanded as in Eq. ([^. Another 
possiblity is to consider non-analytic functions, where this is not possible. The 
simplest of these would be theories with f{R) = . Exact solutions for point- 

masses surrounded by spherically symetric vacua have been found in Refs. 7 & 8. 
A static vacuum solution is given byEI 

= —Ai{r)dt^ + + r^{d9^ + sin^ 6d(j )^), 


Bi(r) 


( 10 ) 


where 




23 ( 1 + 2 ( 5 ) 

V) = T 


.) = 


- cir 


(1-5) 


1 -f Cir’ 


( 1 - 25 + 454 ) 

(1-5) 


(l-25-h4(52)(l-2(5-252) 

A non-static solution, on the other hand, also exists, and is given byll 


where 


ds^ = —A2{r)dt^ + t~ 


/ 1 _ £2 \ 2/o- To \ 4 

Mr) = (j and B^ir) = (l + -) A 

The quantities ci and C 2 that appear in these equations are both arbitrary constants. 
The only other remaining quantitiy is a = 1 — 25 -\- 46 ^. 

Although the weak-field limit of these non-analytic theories is cumbersome to 
probe, these exact solutions can be used to try and construct Swiss cheese models. If 
one attempts this, however, then one finds that no such constructions are possible!^. 
The junction conditions across the boundary cannot be satisfied unless (5 = 0, which 
reduces the theory to that of Einstein. This is true for the geometry described in Eq. 


23(1 + 23 ) 


'B2{r){dr‘^ + r'^{d0^ -\- sin^ OdM) 


( 11 ) 


(T+2<5-1 


(10), and the geometry described in Eq. (11). It is also true of the Schwarzschild 


solution, which is often another solution of the field equations of f{R) theories. 
The only exception to this latter case is if f{R) is linear in R, which is again 
just the Einstein equations with A. So, once again, we find that the only large-scale 
expansion that is possible is the same as that prescribed by the Friedmann solutions 
of Einstein’s equations. 


6. Models of f{R) Dark Energy 

Another family of theories that cannot be expanded as in Eq. ([^ are those con¬ 
sidered in Refs. 9, 10 & 11. These theories have attracted a lot of attention in 
the cosmology community, so it seems well worth investigating them within the 
current context. We start by expanding the quantities / and fa in terms of their 
order-of-smallness in the post-Newtonian expansion. This gives 


/ = M M + 0{e^) and 


/i^ = /r+/f +0(e4), 


( 12 ) 
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where superscripts in brackets denote the order of a quantity in e, which is the 
post-Newtonian expansion parameter. This is exactly the way that the weak-field 
limit of these theories is usually treated. 

If one now tries to build one of the constructions we considered earlier, then one 
is immediately led to the followin^i^ 

/(°)=0 and /n^=qit), (13) 

where q(t) is an, as yet unspecified, function of time only. This leads us to consider 
two possible classes of theory: 


• Class I: These are theories in which q =constant. It immediately follows that 
f{R) — qR + F{R) in this class, where ~ e^. 

• Class II: These are theories in which q ^ constant. This requires q = q{R)^ as 

For this to be true implies i? is a function of time only. 


Let us consider each of these classes separately: 


Class I Theories. This class of theories contains all of those considered in 
Refs. 9, 10 & 11. In this case the field equations 0 and junction conditions 0 
give^ 


^=l(u+ and = ^Fr^ , (14) 

where V^t/ = — 47 rp — F(0)/2 and V^l/ = — 47 rp -|- F(0)/4, and where F{0) = 
constant. The reader can once again note that the U and U functions are exactly 
as one would expect from the solutions of Einstein’s equations with a cosmological 
constant (where F{0) now plays the role of what is usually denoted A). The only 
possible new behaviour comes the the terms containing Fr, which must obej<^ 


V^Fr = ^R- 




(15) 


If, however, we impose the junction conditions from Eq. (|^, then it can be seen that 
the normal derivative of Fr must vanish on the boundary. This means that the extra 
terms in Eq. (14) cannot contribute to the motion of the boundary. We therefore 
find, once again, that the only possible large-scale expansion is indistinguishable 
from that found in the Friedmann solutions of Einstein’s equations with A. 


Class II Theories. The leading-order part of i?, in all theories in this class, 
must be a function of time only. This immediately implies that the leading-order 
part of must also be a function of time only, as = f^'^\R). If this is true, 
it is straightforward to show the post-Newtonian parameter 7 must be given byEH 



1 

2 ■ 


(16) 


This is strongly inconsistent with observational bounds, and occurs because the 
fifth-force that occurs in this class of f{R) theories cannot be screened by the 
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Chameleon mechanism. This can be seen from the equation of motion of the effective 
scalar degree of freedom, which is given by 


= — ^p + homogeneous terms . 


(17) 
( 2 ) 

There does not exist sufficient freedom to create a non-trivial potential for 
in this expression, and so the fifth-force goes un-screened. This means that all 
theories within this class, while potentially admitting interesting new cosmological 
behaviour, cannot be considered observationally viable. 


7. Conclusions 

We find that, in the presence of non-linear structure, the observationally viable 
theories we have studied all evolve on large scale like the Friedmann solutions of 
Einstein’s equations with A. We also note that the oscillatory behaviour that is 
generically present in the Friedmann solutions of many of these f{R) theorie^ 
appears to be supressed. This indicates that either the weak-field limit of these 
theories is destroyed when the oscillations occur, or that the existence of non-linear 
structures prevents the oscillations from occuring at all. Whatismore, in every 
viable case, the effective cosmological constant that emerges must be constructed 
from constant parameters that appear linearly at lowest order in the gravitational 
Lagrangian (he. just like A does). This leads us to the conclusion that these theories 
do not solve any of the problems associated with the cosmological constant, and 
hence that the difficulties created by introducing a new light degree of freedom on 
small scales have no obvious benefit. 
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